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Abstract 

In this paper, we solve a set of hierarchy equations for the re- 
O duced statistical density operator in a grand canonical ensemble for 

an identical many-body fermion system without or with two-body 
interaction. We take the single-particle approximation, and obtain 
^ . an eigen-equation for the single-particle states. For the case of no 

^ | interaction, it is an eigen-equation for the free particles, and solu- 

tions are therefore the plane waves. For the case with two-body in- 
teraction, however, it is an equation which is the extension of usual 
Hartree-Fock equation at zero temperature to the case of any finite 
temperature. The average occupation number for the single-particle 
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states with mean field interaction is also obtained, which has the same 
Fermi-Dirac distribution form as that for the free fermion gas. The 
derivation demonstrates that even for an interacting fermion system, 
only the lowest N orbitals, where N is the number of particles, are 
occupied at zero temperature. In addition, their practical applications 
in such fields as studying the temperature effects on the average struc- 
ture and electronic spectra for macromolecules are discussed. 

PACS 05.30.-d; 31.70.-f; 71.10.-w; 72.80.Le 

1 Introduction 

Quantum statistical mechanics provides a most general theoretical frame- 
work for studying the structure and dynamics of an interacting many-body 
system. Combined with mathematical laws of statistical distribution, it in- 
vestigates the microscopic quantum states of a system with many degrees of 
freedom, its corresponding macroscopic thermodynamic behaviors, and their 
interplay [1-4]. It covers the usual quantum mechanics as a special case at 
zero temperature [5]. 

Experimental investigation of temperature or pressure effects on the elec- 
tronic structure and related spectra for molecules or solids has been for some 
time [6-8]. However, a complete or systematic theoretical work has not been 
developed and seen in the literature. The quantum chemistry at finite tem- 
perature or pressure is not a well-defined or well-established field [9]. In- 
deed, the effect of temperature or pressure on microscopic structure is a 
complicated phenomenon. There exists different functioning mechanisms. 
One consideration is that the variation of the temperature, according to 
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the Fermi-Dirac statistics, will change the thermal probability distribution 
of single-particle states for a free electron gas. Similar situation should be 
expected to occur for an interacting electron system, and therefore its mi- 
croscopic structure will be correspondingly altered. Another consideration 
is that, for molecules or solids, the thermal excitation will cause the change 
of the time scales for the molecular motions. This will most likely bring 
about the transitions of electronic states, and therefore lead to the break- 
down of the Born-Oppenheimer approximation. Electron-phonon interaction 
is a fundamental topic in solid state physics and its temperature dependence 
is well-known. In this article, we tackle the issues in a simpler way. We treat 
only an interacting identical fermion system, or neglect the coupling of the 
electronic motion with those of the nucleus in molecules or solids. We expect 
that some sort of the general conclusions will come out from this study. As 
a matter of fact, this is also the approach usually adopted in non-adiabatic 
molecular dynamics, in which purely solving the eigenequation for the elec- 
trons will provide the reference states for investigating the coupling motions 
between the electrons and the nucleus of the molecules. 

In a recent paper [10], we have deduced a set of hierarchy equations for the 
reduced statistical density operators in both canonical and grand canonical 
ensembles. They provide a law according to which the reduced density matrix 
varies in temperature. They also offer a route for a direct determination of 
the reduced density operator for the statistical ensembles. 

This paper is an initial endeavor, on the basis of the equations we derived 
in paper [10], to investigate the issues related to the interplay between the 
microscopic structure and macroscopic observables, and specifically, the tem- 
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perature effects on the electronic structure of molecules or solids. We take 
the orbital model and solve the equations for the grand canonical ensemble in 
an approximate way. In the next Section, we first consider the case of a free 
electron gas, and evaluate the terms related to the single-particle operator of 
the reduced Bloch equations, we get the solutions of a plane wave as well as 
the usual Fermi-Dirac statistics of a free electron gas. In Section 3, we solve 
the equations for the identical fermion system with two-body interaction. We 
obtain an energy equation for the single-particle states which is the exten- 
sion of usual Hartree-Fock equation at zero temperature to the case of any 
finite temperature. The average occupation number is also obtained which 
has the same form as that for a free fermion gas, which is the Fermi-Dirac 
distribution. The Section 4 is a summary and conclusion. The implications 
and possible applications are discussed. 



2 Free Fermion Gas 



In this and next sections, we consider the solution of reduced Bloch equations 
for the grand canonical ensemble we derived in paper [10] as shown below, 
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and 

h(i) = h(i) - 12. (3) 

The h(i) and g(i,j) are one- and two-body operators, and the rest of sym- 
bols are defined in [10]. We choose a grand canonical ensemble for our study 
because we want to compare our results with those already known for a free 
fermion gas. Just like solving any other equations of motion, the exact solu- 
tions can be obtained only for very few cases. Some sort of approximations 
have to be made. Nevertheless, more universal conclusions might be possi- 
bly obtained. A common approximation scheme is to truncate the hierarchy. 
Depending on the specific physical situation, this is accomplished by approx- 
imating the higher-order density matrix as a functional of the lower-order 
ones. In this paper, the system considered is iV-fermions and the decoupling 
scheme we assume is to express D p+1 by D p in the following way [11-15], 

d p+i = d p AD 1 /D° 

= D 1 A D 1 A ... A D\ /(D°) P , 
p+i 

and 

D p+2 = JJP+ 1 /\ D 1 /D° 

= D 1 A D 1 A ... A D\ / (D°) p+1 . (4) 

p+2 

That is, the pth-order reduced density matrix can be expressed as a p-fold 
Grassmann product of its first-order reduced density matrices. This implies 
that we have taken an independent-particle approximation. To proceed fur- 
ther, we need to evaluate various terms of right hand side of Eq. (1). 
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We first study the case without any interaction, where we only need to 
calculate the first and second terms of Eq. (1). For the first term, 

fiP.D p = H\ ■ D 1 A D 1 A ... A D\ /{D°) p - 1 



v 

I 

i=i 



= AD 1 A... Ah(i) -D 1 A... A D 1 /P°) p_1 - ( 5 ) 



For the second term, the calculation is more complex, and the following 
formulae related to the decomposition of a permutation group in terms of its 
subgroups are needed [16], 

1 



ip+l — 



p+1 
1 

p+l 



A p -J2A p -(p+l,i) 
i=i 
v 



and 



(6) 



(7) 



h(p +l)-A p -(p+l,i)=A p -(p+ 1, i)-]h(i), 

where A p is a pth-order antisymmetric operator, and (p+ 1, i) is an exchange 
between p+l and i. Therefore, 



(p + l)L p p+1 \h(p+l)D p+1 



= (P + l)^ + i 



Z} 1 A. 



p+i 



p+i 



/(D 



0\p 



= (TrhD 1 ) D 1 A D 1 A ... A D\ /(D°) p - 

p 

-- ]T A ... A D l h(i)D l A... AD 1 /(D°) p . 



(8) 



For the special case of p — 1, inserting Eqs. (5) and (8) into Eq. (1) yields 
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This is the Bloch equation for the first-order reduced density matrix of an 
identical particle system with no interaction under single-particle approxi- 
mation. 
Define 

P 1 = D^/D , (10) 
and then Eq. (9) can be simplified as 

-|^ = hp'-p'hp 1 . (11) 

From above equation and its conjugate, we get 

hp 1 -p 1 h = 0. (12) 

This means that h and p 1 commute. Since they are also Hermitian, they 
have common eigenfunctions {\(f>i >}, which are plane waves, 

h\<j>i >= €i\<f>i >, (13) 

and 

p 1 \(f) i >=u([3,p,e i )\<i) l > (14) 

with 

P 1 = I>(/3,A*,ei)|0i >< &|. (15) 

i 

Substituting Eq. (15) into Eq. (11), we can obtain the equation the thermal 
probability u>((3, p, e«) satisfies, 
d 

~ -Qp^iP' ^> € i) = ( £ i ~ P) U (Pi P> e i) - ( e i ~ vWiP, P, ( 16 ) 



Its solution takes the form, 
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which is just the usual Fermi-Dirac distribution. Therefore, by directly solv- 
ing the reduced Bloch equations for the reduced density matrix in a grand 
canonical ensemble under orbital approximation, we not only have obtained 
the eigen solutions of a plane wave but also have recovered the usual Fermi- 
Dirac distribution for the free electron gas. 



3 Hartree-Fock Equation at Finite Tempera- 
ture 

In this section, we solve the reduced Bloch equation (1) for the case of an 
identical fermion system with two-body interaction under the orbital approx- 
imation (4). The first term and the second term remain the same as those for 
the free fermion case. We need to evaluate the remaining two terms. In 
particular, for p = 1, the Eq. (1) reads 

~Y p D l = hD^^^D'-^DW + 2L\ [«K1,2)£> 2 ]+3LS [s(2,3)D 3 . 

(18) 

The last two terms can be evaluated in a straightforward way, which yields 

(19) 



and 



where 



and 



3L l 3 



2L\ [g(l,2)D 2 \ = (J-K)D 1 , 
Tr(gD*) 1 . 



g(2,3)D 



-D\J-K)D\ 



J = Tr 



K = Tro 



9 -D\2;2)\ /D°, 
g-(2,3)-D\2;2)} /D°, 



(20) 
(21) 
(22) 



are called the Coulomb and exchange operators, respectively. The action of 
K is 

K-D\Z;Z) = Tr 2 \g ■ (2, 3) • D l {2; 2)1 / D° ■ D 1 (3; 3) 



Tr 2 



g-D 1 (3;2)-D 1 (2;3)]/D°. (23) 



Substitution of Eqs. (19) and (20) into Eq. (18) yields the Bloch equation 
for the first-order reduced density matrix of N interacting fermions under 
orbital approximation, 

-^ ={F -^ + ^ + T^y-^ {F -^, (24) 

where 

F = h + J-K, (25) 

is called the Fock operator at finite temperature. Redefine the normalized 
first-order reduced density operator 

p 1 = D 1 /D°, (26) 

we can simply above equation into 

-■jlpft = (27) 

From Eq. (27) and its conjugate, we get 

Fp 1 - p 1 F = 0, (28) 

which means that the Fock operator F and p 1 commute. They are also 
Hermitian, and therefore they have common eigenvectors {|0j >}. These 
vectors are determined by the following eigen equation for the Fock operator, 

F\<p t >= e t \<p t > . (29) 



The first-order reduced density operator is correspondingly expressed as 
p 1 = £>(/3,At,e i )|0* >< <t>i\, (30) 

i 

where u(fi, /i, q) is the thermal probability that the orbital is found to be in 
the state {|0j >} at finite temperature T. Substituting Eq. (30) into Eq. 
(27), we can obtain the equation this thermal probability u(/3, /i, ei) satisfies, 

d 

Its solution has the same usual form of Fermi-Dirac statistics for the free 
electron gas as follows, 

u(P,n,€i) = 1 + e ^_ M) , (32) 
with energy levels {e^} determined by Eq. (29). 



4 Discussions and Conclusions 

In this paper, we have solved the set of hierarchy Bloch equations for the 
reduced statistical density operator in a grand canonical ensemble under 
single-orbital approximation for the identical fermion system without or with 
two-body interaction. For the case without any interaction, we not only get 
a plane wave solution for the eigenstates and eigenvalues, but also recover 
the usual Fermi-Dirac distribution of a free electron gas. For the situation 
with two-body interaction, we obtain an eigen-equation for the single-particle 
states. It is the extension of usual commonly used Hartree-Fock equation at 
the absolute zero temperature to the case of any finite temperature. The 
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average occupation number formula for each single-particle state is also ob- 
tained, which has the same analytical form as that for the free electron gas 
with the single-particle state energy determined by the Hartree-Fock equa- 
tion at finite temperature (29). 

From Eqs. (21), (22) and (30), we see that the Coulomb operator J, the 
exchange operator K, and therefore the Fock operator F are both the co- 
herent and the incoherent superpositions of single-particle states. They are 
all temperature-dependent through an incoherent superposition factor, the 
Fermi-Dirac distribution, u>(/3, n,€i). Therefore, the mean force and corre- 
sponding microscopic structure are temperature-dependent. 

It is easy to see that, when temperature is zero, only N single-particle 
states with energy levels small than the chemical potential /i are occupied 
as for the free electron gas. That is, the fact that only lowest N orbitals 
(or holes) are occupied holds also for the interacting case. Therefore, only 
a single determinant wavefunction is enough for a description of interacting 
iV-particle wavef unctions. Say in another way, there exists a correspondingly 
effective and strict one-particle state description of an interacting iV-particle 
system. This is also the spirit of current density functional theory [17], and 
is consistent with the third law of thermodynamics. 

Even though we have simplified the structure issue for a molecule or a 
solid at finite temperature, the equations we obtained in this paper will find 
wide applications. The examples include the investigation of temperature 
effects on the average structure and the electronic spectra of macromolecules 
and the study of conduction electrons in metals and so forth [18-21]. 
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